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Abstract 

In earlier work [2^ the authors obtained integral formulas for probabilities 
for a single particle in the asymmetric simple exclusion process. Here formulas 
are obtained for joint probabilities for several particles. In the case of a single 
particle the derivation here is simpler than the one in the earlier work for one 
of its main results. 

I. Introduction 

The asymmetric simple exclusion process (ASEP) is a continuous time Markov 
process of interacting particles on the integer lattice Z. Each particle waits expo- 
nential time; then with probability p it jumps one step to the right if the site is 
unoccupied and otherwise it stays put; with probability g = 1 — p it jumps one step 
to the left if the site is unoccupied and otherwise it stays put. We refer the reader to 
Liggett [m |T2] for a precise definition of the model. 

When particles can jump only to the right or only to the left, the model is referred 
to as the totally asymmetric simple exclusion process (TASEP). This process is quite 
special in that it, unlike ASEP, is a determinantal process [H [131 [221 [23] . Fo^' example, 
in A^-particle TASEP the transition probability from an initial state Y = {yi, . . . , y^} 
to a state X = {xi, . . . ,X]^} in time t is given by an x determinant [22]. For 
TASEP with step initial condition {yi = i) Johansson, in the seminal paper [13], 
proved a limit law for the particle current fluctuationsll] 



^The limiting distribution is the distribution function F2 of random matrix theory. 
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Subsequently, these particle current limit laws for TASEP were extended to other 
initial conditions including periodic initial conditions [5], stationary initial conditions 
|9] and the general two-sided Bernoulli initial condition [Sj |9] . Much recent work on 
TASEP has focused on the joint distributions of the particle current [21 El 13 El [HI 
HHIIIT]. In particular for step initial condition, the joint distribution of the associated 
height function and its KPZ scaling to the Airyg process [lll[T8] is now well developed 
for TASEP [15] . We refer the reader to [10] for a recent review of TASEP and closely 
related random growth models. 

For ASEP the results have been restricted to probabilities for a single particle 
[211 [251 (211 [23 • III [23 [21] the present authors extended Johansson's limit law for 
TASEP to ASEP thus proving a stronger form of KPZ universality. Using results 
from [25], Sasamoto and Spohn [191 [201 [21] and Amir, Corwin and Quastel [1] have 
obtained explicit formulas for the exact height distribution for the KPZ equation 
with narrow wedge initial condition. This height distribution interpolates between a 
standard Gaussian distribution for small time and the F2 distribution for large time. 
To extend the results of [H [191 [201 121] and show that the Airyg process is the large 
time limit of KPZ with narrow wedge initial condition we would first try to extend 
the results of [211 ES] to joint distributions of particle positions in ASEP. In this 
paper we take a first step in this direction in that we derive integral formulas for joint 
probabilities. Whether these formulas can be used to study the KPZ scaling limit 
remains to be seen. 

Here is an outline of the paper. We consider an A^-particle system with initial con- 
figuration Y = {i/i, . . . , Un}, where yi < ■ ■ ■ < y^, and denote by {xi{t), . . . , xp^(t)} 
the configuration of the system at time t. Our starting point is Theorem 2.1 of |24j . 
a formula for the probability 

FY{xi{t)=Xi, i = l,...,N) (1) 

(the subscript indicating the initial configuration), which is a sum of A^! multiple 
integrals, one for each permutation in the symmetric group S^y. The integrals are 
over small circles about the origin. In Sec. II we obtain a formula for 

Py(xi(t) =Xi, ^ = l,...,m), (2) 

the joint probability for the first m particles, by summing ([1]) over all Xi with i > m 
which satisfy Xm < Xm+i < ■ ■ ■ < xn < 00. We may do this because the integrals are 
over small circles. Application of combinatorial identity (1.6) of [21] replaces the sum 
of A^! integrals by a sum of N\/{N — m + 1)! integrals. This is Theorem 1. 

To obtain from this a formula for arbitrary consecutive particles, 

Fyixiit) = Xi, i = n, . . . ,m), (3) 
^KPZ refers to Kardar, Parisi, and Zhang [T5] . 
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we have to sum over all Xi with i < n which satisfy — oo < xi < • • ■ < Xn-i < and 
to do this requires first replacing integrals over small circles by integrals over large 
ones. This is accomplished with the help of Lemma 3.1 of |24], and gives Theorem 2. 
Then summing over the indicated Xi and applying combinatorial identity (1.7) of [21] 
give a formula for ([3]) as a sum of integrals over large circles. This is Theorem 3. 

One could obtain a formula for the general joint probability 

by taking n = mi and m = rrir in Theorem 3 and summing over those Xj satisfying 
Xm < Xi < Xm +1 • These are finite sums and lead to a rather complicated expression 
which we shall not write down. There does not seem to be a combinatorial identity 
that simplifies it. 

Finally we consider the special case Py(xm(t) = Xm), the probability for a single 
particle. The formula for consecutive particles simplifies by use of combinatorial 
identity (1.9) of [23], and yields Theorem 4. This is exactly Theorem 5.2 of [23]- 
The proof here is much simpler than the one in [23], although they use the same 
ingredientsif] The awkwardness of that proof made it seem unlikely at first that one 
could find formulas for joint probabilities without great effort. 

Although we assumed that the initial configuration Y was finite. Theorems 2-4 
extend to configurations that are semi- infinite on the right, in particular ioY = Z"*" 
(step initial condition). This is shown in the last section. 



II. The first m particles — small contours 



Theorem 2.1 of [24J is the formula, valid when p 7^ 0, 

Py(Xi(t) =Xi, ^= 1,...,A^) 

E / • ■ 7 n c» n ^ '^'-^ 'di,... di^Si (4) 



where 



A^ = sgna ^.^ . /(e,n=P + ?a'-e, siO=Pr' + q^-l- (5) 

i<j 



•^The theorem was the starting point for some of the advances mentioned above, and so it is useful 
to have this more straightforward proof. 

''in general the product symbol Y[ refers only to the first term after it. 



3 



Here Cr is a circle with center zero and radius r, which is so small that all nonzero 
poles of the integrand lie outside Cr- (All contour integrals are to be given a factor 
l/27ri) 

To obtain Py(xi(t) = Xj, i = 1, . . . , m) we sum over all Xm+i, ■ ■ ■ ,xn such that 
m+1 < ■ ■ ■ < xat < cxD, which we may do when r < 1. The result of the 
summation is that the product ^^^^^^ in (jlj) is replaced by 

<icj(m+l) SCT(m+2) ■ ■ ■ ?(T(Af) 



,.,,..-rv.-{-+2)---^<x{iV))---(l-e<x{JV))' 

We take a fixed ordered (m — l)-tuple A = . . . ,irn-i} with distinct ii < N 
and consider first the sum over all permutations a such that 

o-(l) = £i, . . . , a(m - 1) = im-i- 
(This meaning of A as an (m — l)-tuple will be retained throughout.) 
The product U^<^^Q^) U^>m Cil) becomes 

(in the second factor here, and elsewhere below, A is thought of as a set) and the 
numerator in the expression for Acr may be written 

j¥=h j¥=h:i2 m<i<j 

iGAjjr^A i<j m<i<j 

As for sgn cr we will consider a new guise for a as its restriction to A'^ = [1, iV]\A. 
(It may be associated in an obvious way with a permutation in STv-m+i, as A may be 
associated in an obvious way with a permutation in S^-i-) What we now call sgn cr 
will be different from the old sgn cr, because the number of inversions in the original 
cr equals the number of inversions of the new guise a plus the number of inversions 
in the original cr which involve some ii. This number is 

"y^^jii — i) + # of inversions in A. 

Thus the original sgn cr in ([5]) is to be replaced by 

(_l)Ei<m(^i-0 sgn A sgn cr. 
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where now a is in its new guise. 

Identity (1.6) of ^24] tells us that, with the new guise a, 



J^sgna Y[ /(^^(.), ^^(^.)) —_- ^ 



/• ^-2 (:N—m 

?o-(m+l) So-(m+2) ' ' ' ^o-(Af) 



n - 



^{N-m+l){N-m)/2 /' ]^ _ JJ" ^ A 



11(1 



Now we sum over all ordered (m — l)-tuples A and obtain 
Theorem 1. We have when p 

Py(x,(t) = x„ z = 1, . . . , m) = ^(-1)J:.<-(^»-*) sgn a p(A^— +i)(a^— )/2 

A 



n - «•) 

i^A 11^ ^■J'' i^A i 



III. The first m particles — large contours 



Here we shall find another representation of the same probability in which all 
contours of integration are Cr, where R is arbitrarily large. Observe that some factors 
/(^j, ^j) in the denominator in the integrand are cancelled by factors in the numerator. 
We assume at first that p, g 7^ 0. 

Step 1. We show first that we may take all ^^-contours with £ G A to be Cr. Let 
us first expand the ^^^-contour, where £a = max A. There appear to be poles at 



1 - 

^From now until the last section N — \Y\, all indices belong to [1, N], and all sets of indices are 
subsets of [1, N]. 
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coming from the factor /(^^j,, ^j) in the denominator. Since ^a = max A we must have 
j ^ A. But then this factor is cancelled by the same factor in the numerator. So the 
only poles when we expand the ^^^-contour are at 

which comes from the factor /(^j, in the denominator when i < i^- 

We show that the integral with respect to C,i of the residue at this pole equals 
zero. When we make the substitution — )■ (^j — p)/q^i we find that e{C,i) + e{^e^) 
becomes analytic at C,i = 0. As for the first quotient in the integrand, each factor / 
having as one of its variables is of the order as — )■ 0. There are — 1 such 
factors in the numerator and A^ — 2 in the denominator (since we don't include the 
factor /(^j, ^i^) that gave rise to the pole). Thus this quotient is 0(^~^). The residue 
of at the pole is so these combine to give the power ^j^'^. As for 

the product of powers of the variables at the end of the integrand, when i ^ A the 
product of those involving C,i is 

since Xa < Xm, and i < ia- Thus the integrand is analytic inside Cr and so its integral 
is zero. If i G A then i = with (3 > a. (Otherwise the factor /{^i^, ^e^) in the 
denominator is cancelled by the same factor in the numerator.) Then the exponent 
Xm above is replaced hj x^, and Xj3 > Xa- 

We have shown that we may expand the ^^^-contour to Cr. Next we expand the 
^^^-contour where is the second-largest element of A. The only difference from 
what went before is that there could be a pole at 

coming from the factor /(^£^, $,£^) in the denominator. But G Cji, and if R is large 
enough the pole would lie inside Cr and so would not be passed in the expansion of 
the ^^^-contour. Thus we may expand the ^^^-contour to Cr 

Continuing in this way we see that all the ^rcontours with £ G A may be taken 
to be Cr. 

Step 2. Next we expand the ^j-contours with i ^ A, and when we do this we 
encounter poles at = 1 as well as those coming from the first quotient in the 
integrand. As above the latter poles will not contribute, but those from = 1 will. 
An application of Lemma 3.1 of [23] will tell us the result. 
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We restate the lemma to make it compatible with the present notation. Let T be 
a finite set of indices and let 

9 = 9{ii)i&T 

be a function that is analytic for all 7^ 0. Assume that for i > k 
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as ^fc — )■ 0, uniformly when all the with j ^ i, k are bounded and bounded away 
from zero. For S G T denote hj gs the function obtained from g by setting all C,i 
with i ^ S equal to 1. (In particular g = gjn.) Define 



T ((:) ^ TJ ~ 9s{0 



Then when p, q ^ 0, 

where is so large that all the zeros of the denominators lie inside Cr. (When S is 
empty the integral on the right side is interpreted as (/(I, . . . , 1).) Here t = p/q and 

a{S, T) = #{(*, JET, t>j}. 

We apply this with T = and 

9{^^)^^A^ = (1 - n ^0 n e^-A-(«^)* 



j<m ieA ieA 



z or J £ A 



The integral on the left side of ([7j) is then the integral in Theorem 1. That g is 
analytic for ^ despite the denominator in the integrand follows from the fact 
that R may be arbitrarily large. The argument leading to (|6]) is like the analogous 
argument in Step 1. There are m — 1 factors of order in both the numerator and 
denominator of the quotient in the integrand, and the product ik^''~^ becomes 
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since i > k. 

So the lemma applies. If \S\ = k, then the coefficient on the right side of is 

^_(N-m+l)(N-m)/2 ^a{S,K'')-k ^k{k-l)/2_ 

However there is another factor arising from 

n /te'^^) 



n /te'^.) 



i or J £ A 



when any with k G A'^\5' is set equal to one. This factor is 

- i'_^^#{ieA:j>fc} 



Hence the total factor when the with k G A'^\S are set equal to one is 

If we use the bilinearity of (j{U, V) in its two variables, and the easy fact 

aiU, V) + aiV, U) = \U\ \V\ + \UnV\, 

we see that 

a{A, S"" n A") = a{A, A") - a{A, S) = J2^^- - l)/2 + (^{S, A) - (m - l)k. 
Thus (IE]) times ([9]) equals 

where we have written (t{S) for a{S, A) + cr{S, A'^) = (t{S, [1, A^]). 

We replace each integral in the formula of Theorem 1 by the sum over 5 C A'^ we 
obtain using ([7]). The result is 



8 



Theorem 2. We have when g 7^ 0, 

A 
SCAC 

n - 

jes 11^ ^i'' «<"^ jes i i 

The summation is over all sets S and (m — l)-tuples A disjoint from S. Here = IS*! 
and in the integrand indices that are not specified belong to 5* U A. 

We assumed at first that p, g 7^ since that was required by Lemma 3.1 of [21]. 
To obtain the relation when p = by passing to the p — )■ limit we need only 
observe that the power of r in the coefficient is nonnegative and that the integrand 
is continuous at p = when R> 1. 



IV. Consecutive particles 



Here we find a formula for Py(a;j(t) = x,, i = n, . . . , m). If we take the preceding 
formula and sum over all xi, . . . , with —00 < Xi < ■ ■ ■ < x„_i < x„, which we 
may do when i? > 1, the partial product H^n'^r/ replaced by 

(61 • ■ ■^^n-l)"''" 



Our (m — l)-tuple A may be written in the obvious way as {Ai, A2}, where 

Ai = {£1, . . . , C-l}, A2 = {C, . . . , f-m-l}. (11) 



Until now we did not differentiate in the notation between an ordered set and the 
corresponding unordered set. Which was meant was always clear from the context. 
Now we introduce a notation when we want to distinguish the two: for an ordered 
set A we denote by X(A) the corresponding unordered set. 
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In the integrand in Theorem 2 the only factor that depends on A and not just on 
the set X(A) is 

i<j i<j<n i<'n',j'>n n^i<j 

The first factor depends on Ai, the second depends only on X(Ai) and X(A2), while 
the third depends on A2. 

In Theorem 2 we are to sum over all sets S and all (m — l)-tuples A disjoint from 
S. Let us fix S first and sum over all A disjoint from S. This sum over A is the 
sum over all Ai and A2 as in ( ITTi) . We do this sum by first taking a fixed set Si 
disjoint from 5* with 15*11 = n — 1, and summing over all Ai with X(Ai) = Si. Then 
we must sum over Si and, of course, A2 (which is disjoint from and S) and S. To 
recapitulate: the summation in the statement of Theorem 2 may be replaced by 

E 

S, S'i,Ai,A2 

where S and Si run over all disjoint sets with IS*!] = n — 1, and Ai and A2 run over 
all (n — l)-tuples and (m — n)-tuples respectively satisfying 

X{Ai) = Si, X{A2)c{SuSiy. 
With this notation the coefficient in Theorem 2 may be written 

Notice that Ai appears here only in the factor sgnAi. The sum of those terms 
involving Ai as distinguished from 5*1 is 

(Recall the factor (fTOj) .) Identity (1.7) of [2l] tells us that this sum equals (^^^ ■ ■ ■ J^" 
times 

n (^^- - n (^^- - 

^(n-l)(n-2)/2 '-jegi ^ g("-l){"-2)/2 '-jegi 



jeSi jesi 

We write the last product here times the analogous product in the integrand in 
Theorem 2, where S appears instead of 5*1, as 



n -^^ 



,.:<^,o. (_l)-(5i,5) 



i, jesuSi 



n (1-^^) n (^.-^^ 

j&suSi ieSi, jes 
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Putting all these things together and using a{S, Si) + (t{Si, S) — {m — l)k we 
obtain 

Theorem 3. We have when g 7^ 0, 

FY{xi{t) = Xi, i = n,. . .,m) 

^ ^ (^_-^y{SuSi,A2)+n-l sgnTV^ -r^{'S'U5iUA2)-m(m-l)/2-mfe ^fc(A;-l)/2+(n-l)(n-2)/2 
S,5i,A2 

n fi^i:Q n n /(^.cj 

^ / / ieSi,jeSuA2 i£A2,jeS n<i<j 

n Q 

i<j 



n (^^■-^^) 



i,jesuSi 



^ n (1-^^) n (^^•-^') 

iSSi n<i<m j65 i i 

The summation is over all disjoint sets S and Si with l^*! | — n — 1 and (m — n)-tuples 
A2 disjoint from S U S*!. Here k — \S\ and in the integrand indices that are not 
specified belong to 5" U 5*1 U A2. 

V. One particle 

In this case n = m and so A2 disappears, the first line of the right side becomes 

^<T(SuSi)-m(m-l)/2-mA: ^fc(fc-l)/2+(m-l)(m-2)/2 

S,Si 

and the factor 

n^ " n c n^r 

ieSi n<i<m jeS 

becomes 

n 

ieSuSi 

Now we take a fixed set 5*3 with l^al = m + k — 1 and sum over all partitions 
Ss = {S, Si} with \S\ = k, |5'i|=m — 1. The only part of the sum that depends on 
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S and 5*1 individually is 



n 



S, Si 



i€Si, jes 



If we observe that S is the complement of 5*1 in 5*3 we see that we can apply identity 
(1.9) of [21] (with S there 5*1 here, with m there m — 1 here, and with there m+k — 1 
here), which tells us that the sum equals 



(m-l)fc 



m + k — 2 
m—1 



where the r- binomial coefficient [^]^ is defined by 

(1 - r^) (1 - r^-i) ■ ■ • (1 - r^-+i) 



n 



'l-r)fl-r2 



'1 - T^- 



Hence, after some algebra, our formula becomes 
Theorem 4. We have when g 7^ 0, 



n 



t53|>m 



m — 1 



where all indices in the integrand run over S3. 
This is exactly Theorem 5.2 of |24j . 



VI. Semi-infinite configurations 



Here we show that Theorem 2, and consequently also Theorems 3 and 4, hold 
when Y is semi-infinite on the right, if we take the sums over finite sets 5* C A^. It 
follows from the fact that ASEP is a Feller process [11] that the probability for the 
semi-infinite initial configuration Y equals the N ^ oc limit of the probability for 
initial configuration {yi, . . . , yjv}. Thus we sum over only those Ac [1, A^] and those 
S" C [1, A^] in Theorem 2 and then pass to the limit. The limit will be the sum over 
all A and (finite) S if the resulting series is absolutely convergent. This we now show, 
using the fact that R may be taken arbitrarily large. 
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Consider the integrand first. Each / factor is of the order B? for large R while 
each is of order R. Combining the estimates for all factors in the integrand other 
than Wi^i'^"'^ gives, after a little algebra, 0{R~^^/'^^^^^^). (Recall that m is fixed.) 
Since each ?/i > + i — 1, we have Y\i — 0{R~'^^^^^^~^'^^'^^) . Thus the integrand 



li Sj 

is 

-fc2/2-a(AU5)+0(fc) 



0{R' 

Another factor R^ comes from the domain of integration, but this does not change 
the bound. 

The coefficient on the right side of Theorem 2 has absolute value at most r°'i^'~>^)+oik) _ 
If we take i? > r^, as we may, then this combined with the preceding bound gives 

This a bound for the summand in Theorem 2. Now we show that the sum of these 
over all A and S is finite. For any s we have 

#{(A, S) : a{A U S) = s} < 2' s'^'K 

The reason is that if cr(A U 5*) = s then the largest element in A U S* is at most s. 
So there are at most 2^* choices for S, snd having chosen S there are at most s*""^ 
choices for A. It follows that the sum of f|T2l) over all (A, S) is at most a constant 
times 

2S g-m-l ^-fc2/2-s/2+0(fc) 

fc,s>0 

which is finite when i? > 4. 
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